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ABSTRACT 

This  report  summarizes  my  research  on  excitation  and  propagation  of  regional 
waves,  mainly  Pn,  in  the  past  year.  The  nonlinear  spectral  inversion  method  of  Xie  (1993) 
was  modified  to  better  estimate  source  and  path  spectral  parameters  using  Pn  and  Lg. 
The  modification  was  made  to  further  reduce  effects  of  deviations  of  sources,  paths  and 
sites  from  idealized  theoretical  models  used  in  the  inversion.  An  extensive  attempt  was 
also  made  to  estimate  the  Pn  geometrical  spreading  term  in  reasonable  ID  structures  for 
central  Asia.  Synthetic  calculations  of  Pn  propagation  in  central  Asia  were  conducted 
using  three  recently  published  ID  velocity  models.  The  results  typically  show  that  a 
A^^Aq/A”,  with  m  ranging  between  1.3  to  1.5,  is  roughly  valid,  but  one  should  expect 
localized  deviation  of  Pn  amplitude  from  that  model  by  at  least  a  factor  of  2. 

We  have  collected  more  than  one  hundred  Pn  spectra  from  broad-band  seismic  sta¬ 
tions  in  central  Asia.  Pn  amplitude  across  the  Kyrghistan  network  (KNET)  from  sources 
located  in  and  around  the  Lop  Nor  Test  Site  is  found  to  be  highly  unstable.  Despite  the 
small  aperture  of  the  KNET  (~  2  °),  the  time  domain  Pn  amplitude  varies  by  as  much  as  a 
factor  of  19.  Spectral  analysis  reveals  that  the  amplitude  variation  primarily  occurs  at  low 
frequencies  (~  Ihz).  Due  to  the  pronounced  topographic  variation  near  the  stations,  fre¬ 
quency-dependent  polarization  analysis  yields  little  information  on  the  origin  of  the 
amplitude  variation.  Analysis  of  Pn  propagating  through  the  network  in  different  direc¬ 
tions,  on  the  other  hand,  suggests  that  the  amplitude  variation  is  not  due  to  a  systematic 
variation  in  site  responses.  The  more  likely  cause  of  the  amplitude  variations  is  deep- 
seated  3D  structural  complications,  such  as  a  topographic  variation  at  the  Moho.  We 
demonstrate  that  a  2D  Moho  topography,  likely  associated  with  the  drastic  tomographic 
variation  in  the  earth’s  surface  east  of  the  KNET,  could  have  preduced  a  ray  pattern  that 
quantitatively  gives  rise  to  the  observed  Pn  amplitude  anomaly. 

Pn  amplitudes  at  KNET  stations  from  some  Northern  Xinjiang  earthquakes  are  much 
weaker  than  from  the  nearby  Lop  Nor  explosions,  making  Pn/Lg  ratio  an  effective  dis¬ 
criminant  between  these  particular  earthquakes  and  underground  nuclear  explosions.  The 
weak  Pn  amplitude  from  these  earthquakes  are  likely  caused  by  a  null  in  the  Pn  radiation 
pattern  of  the  earthquakes  involved,  whose  mechanisms  are  likely  dominated  by  an  east- 
west  strike  slip  motion.  Therefore,  it  might  be  misleading  to  extrapolate  the  low  Pn/Lg 
amplitude  ratio  to  earthquakes  with  other  source-station  configurations.  For  the  earth¬ 
quakes  that  occurred  in  southern  Tien  Shan,  for  example,  observed  Pn  amplitudes  at  the 
KNET  stations  are  typically  quite  high.  This  is  consistent  with  higher  Pn  radiation  pat¬ 
terns  predicted  for  focal  mechanisms  that  are  predominantly  reverse  slip,  which  is  typical 
for  earthquakes  in  southern  Tien  Shan. 

Research  so  far  demonstrates  that  the  Pn  amplitude  in  central  Asia  are  affected  sig¬ 
nificantly  by  3D  structural  complications  and  source  radiation  patterns.  Any  attempt  to 
utilize  Pn  in  discriminating  explosions  from  earthquakes,  or  in  characterizing  source/path 
properties,  must  account  for  these  complications.  An  extensive  analysis  is  underway  to 
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compare  the  stability  of  spectral  inversion  using  Pn  alone  with  that  jointly  using  Pn  and 
Pn  coda  (or  Pg),  which  is  expected  to  be  more  stable  than  Pn. 
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1.  INTRODUCTION 

Lg  and  Pn  are  two  of  the  most  prominent  regional  phases  used  in  the  discrimination 
of  explosions  from  earthquakes.  Xie  (1993)  developed  a  nonlinear  inverse  method  to 
simultaneously  estimate  the  seismic  moment  (Mq),  comer  frequency  (fc)  and  path- 
variable  Qo  and  77  (Q  at  1  hz  and  its  power-law  frequency  dependence,  respectively)  using 
the  Lg  spectra  from  both  explosions  and  earthquakes.  The  method  was  applied  to  21 
underground  nuclear  explosions  (Xie,  1993;  Xie  et  ai,  1996)  and  53  earthquakes  (Cong 
et  ai,  1996)  in  central  Asia.  Major  findings  of  these  studies  include  (a)  The  Lg  Q  values 
at  1  Hz  for  numerous  paths  in  central  Eurasia  were  found  to  agree  rather  well  with  those 
predicted  using  a  tomographic  Lg  coda  Q  map  (Xie  &  Mitchell,  1991).  (b)  For  both 
underground  nuclear  explosions  and  the  earthquakes  studied,  the  logarithm  of  Lg  Mq  val¬ 
ues  correlate  linearly  with  the  body  wave  magnitude  (M^),  with  slopes  of  slightly  greater 
than  1.0.  (c)  The  Lg  Mq  values  tends  to  scale  scale  with  with  the  value  of  a  ranging 
between  about  3.6  to  4.0.  (d)  The  scaling  between  Lg  Mq  and  fc  values  for  the  explo¬ 
sions,  estimated  using  the  earthquake  source  model  (i.e.,  the  o)^  source  model  without 
spectral  overshoot,  often  known  as  the  Bmne’s  model),  differ  from  that  between  the  Lg 
Mq  and  values  for  the  earthquakes.  The  main  difference  between  the  two  scalings  is 
that  for  a  given  Mq,  the  explosions  tend  to  have  higher  values.  This  suggests  that  the 
Lg  from  explosions  has  a  richer  high-frequency  content,  as  compared  to  Lg  from  earth¬ 
quakes.  This  phenomenon  appears  to  be  systematic  in  central  Eurasia,  and  appears  to  be 
opposite  to  that  observed  in  the  western  U.S.  It  thus  raises  a  question  on  the  transportabil¬ 
ity  of  any  explosion  discriminant  based  on  the  Lg  spectral  content  alone. 

The  main  objective  of  this  research  is  to  extend  the  inverse  method  of  Xie  (1993)  to 
analyze  Pn  spectra  fi'om  both  explosions  and  earthquakes.  The  analysis  should  yield  Pn 
source  Mq,  fc,  and  path  Qq,  77  values.  These  values  will  advance  our  fundamental  under¬ 
standing  on  how  the  excitation  and  propagation  of  Pn  differ  from  those  of  Lg.  By  com¬ 
paring  the  Mq  ~  fc  scalings  derived  using  Pn  from  explosions  and  earthquakes  with  those 
derived  using  Lg,  we  will  investigate  if,  when  and  why  the  Pn/Lg  ratio  works  as  an  dis¬ 
criminant  between  explosions  and  earthquakes. 

This  research  is  composed  of  several  tasks.  The  first  is  to  develop  a  relatively  reli¬ 
able  and  robust  method  to  measure  the  source  Mq,  fc,  and  path  Qq,  tj  values  using  the  Pn 
phase.  This  task  is  difficult  not  only  because  of  the  complication  in  estimating  Mq,  fc  val¬ 
ues  due  to  the  source  model  dependence,  as  we  have  discovered  for  the  inversion  of  the 
Lg  spectra,  but  also  because  of  the  effects  of  non-isotropic  source  radiation  pattern  and 
3D  path  structural  complexities  in  the  Pn  spectra.  The  second  task  of  this  research  is  to 
analyze,  with  the  method  developed,  Pn  spectra  from  many  events  to  obtain  their  Mq  and 
fc  values,  as  well  as  the  Pn  Qq  and  77  values  for  various  regions  such  as  central  Asia  and 
middle-east.  The  third  task  is  to  compare  the  Mq  ~  fc  scalings  for  earthquakes  and 
explosions,  derived  using  Pn,  with  those  derived  using  Lg,  to  obtain  the  differences  in 
those  scalings.  In  this  report,  we  summarize  major  accomplishments  of  this  research 


1 


achieved  over  the  past  year. 


2.  MODIFICATION  OF  THE  INVERSE  METHOD 


§2.7  Practical  problems  affecting  the  spectral  inversion 

High-frequency  regional  seismograms  are  affected  by  many  factors,  ranging  from 
non-isotropic  source  radiation  pattern,  3D  structural  complications  along  the  paths,  site 
amplification  and  occasional  instrument  miscalibrations.  Quantifying  these  effects  in 
deterministic  manner  is  impractical  due  to  the  many  degrees  of  freedom  involved.  To 
study  Lg  spectral  characteristics,  Xie  (1993)  adopted  the  stochastic  modeling  used  by 
Street  et  al  (1975),  Herrmann  and  Kijko  n  (1983),  Xie  &  Nuttli  (1988)  and  Sereno  et  al. 
(1988).  In  the  stochastic  modeling  A,(/),  the  Lg  or  Pn  spectra  at  an  /th  station,  is  mod¬ 
eled  as 


A-iif)  =  S(f,  0i)G(Ai)  exp 


r 


VMi(f) 


W) 


(1) 


where  /  is  frequency,  9i,  A„  and  Qi(f)  are  the  azimuth,  distance  and  quality  factor  of  Lg 
or  Pn  from  the  source  to  the  ith  station.  Vg  is  the  mean  group  velocity  of  Lg  or  Pn,  X,  (/) 
is  the  site  response  at  the  ith  station,  and  S(f,  ^,)  is  the  source-to-Lg  (or  Pn)  spectral  radi¬ 
ation  in  the  direction  of  or  source  spectrum,  which  can  be  represented  by 


MoR(0i)  1 

4^  pv^  1  +  ff/fc 

MoRiBi)  _ 1 _ 

[l  +  (1  -  2ll)fV  f}  +  yS2/V  /< 


for  earthquakes 
for  explosions 


(2) 


where  Mq  and  are  seismic  moment  and  comer  frequency.  For  Lg  spectra  p  and  v  are 
the  averaged  cmstal  density  and  shear  wave  velocity;  for  Pn  they  are  the  source-zone 
density  and  compressional  wave  velocity,  respectively.  R{9i)  is  the  azimuthally  varying 
source  radiation  pattern.  In  the  above  equation,  we  have  used  the  omega-square  model 
without  overshoot  for  the  earthquakes  (often  referred  to  as  the  "Bmne’s  model")  and  the 
omega-square  model  with  overshoot  for  the  explosions  (the  modified  Mueller-Murphy 
model;  Sereno  et  al.,  1988).  G(A,)  in  Eq  (1)  is  the  geometrical  spreading  factor  and  typi¬ 
cally  has  the  form 

G(A,)  =  A;'(VAi)”  (3) 

where  Aq  is  a  reference  distance,  and  m  is  the  decay  rate  of  Aj(/)  at  large  distances 
(A  >  Aq).  Typical  values  of  Aq  and  m  are  about  100  km  and  0.5  for  Lg  (Street  et  al, 
1975),  and  about  1  km  and  1.3  for  Pn  (Sereno  et  al,  1988).  In  the  next  section  we  will 
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further  discuss  the  values  of  Ao  and  m  for  Pn  in  central  Asia.  0,(/)  in  Eq  (1)  is  the  appar¬ 
ent  Q  (quality  factor), 

QiU) = Qoi  r‘  (4) 

where  Qo/  and  rji  are  apparent  Q  to  the  ith  station  at  1  Hz  and  its  power-law  frequency 
dependence,  respectively.  To  develop  a  practical  algorithm  for  inverting  spectral  charac¬ 
teristics  of  regional  wave  excitation  and  propagation,  Xie  (1993)  made  an  important 
assumption  that  for  a  given  explosion  event,  R{6i),  Ao  and  m  in  Eq  (2)  and  (3)  are  not,  but 
Qqi  and  77,  are,  station-variable.  Cong  et  al.  (1996)  extended  this  assumption  to  Lg  from 
earthquakes.  When  estimating  Mq,  fc,  Qoi  and  tji  values,  this  assumption  may  well  cause 
bias.  However,  I  think  that  for  modem,  broad-band  seismic  stations,  a  reasonable  estimate 
of  the  variation  in  Lg  spectra  caused  by  the  combined  effect  of  the  station  variable 
X(f),  Aq  and  m  and  the  instmment  malfunction  is  within  about  a  factor  of  2.  At  large  epi- 
central  distances  (roughly  A  >  2.8  Vg  Q{f)  /  {n /)),  the  error  caused  by  this  effect  is  vir¬ 
tually  ignorable.  At  smaller  distances,  however,  there  are  indications  that  this  effect  may 
cause  significant  error  in  Qq  and  tj  estimate.  Cong  et  al.  (1996),  for  example,  did  not  use 
Lg  spectra  from  stations  with  smaller  distances  (<  about  650  km)  since  the  resulting  Qq 
and  1]  values  are  too  inconsistent  with  those  previously  found  for  central  Asia. 

Another  source  of  error  of  the  spectral  inversion  of  Lg,  Pn  and  any  other  high- 
frequency  phases  comes  from  the  tendency  for  individual  source  spectra  to  deviate  from 
the  idealized  stochastic  models  in  Eq  (2).  Xie  et  al.  (1996)  found  that  the  values  of  Mo,  fc 
resulting  from  inverting  Lg  spectra  of  21  underground  nuclear  explosions  are  dependent 
on  which  of  the  source  models  in  Eq  (2)  was  used,  despite  that  the  high-frequency 
asymptotes  of  the  two  models  were  similar.  Recently  Xie  (1996)  also  found  that  the  Mq, 
fc  values  estimated  for  the  1995,  Western  Texas  earthquake  (M^y  ~  5.7)  are  dependent  on 
the  source  model  used.  Since  the  estimates  of  Mq,  fc  and  path  Q  values  are  dependent  on 
the  source  model  used,  a  deviation  of  the  source  spectra  of  any  event  from  the  idealized 
source  model  used  in  a  spectral  inversion  will  undoubtly  cause  bias  in  the  estimated  Mq, 
fc,  Qq  and  77  values. 

These  biases  are  expected  to  be  more  serious  when  one  deals  with  Pn.  This  is 
because  that  (1)  The  distance  range  for  Pn  observation  is  between  about  2  to  12°,  which 
is  significantly  narrower  than  about  2  to  35°,  the  distance  range  for  Lg  observation. 
Therefore  the  condition  of  A  >  2.8  Vg  Q(f)  /  (trf)  is  less  frequently  satisfied  for  Pn.  (2) 
Pn  is  primarily  comprised  of  only  one  ray,  as  against  many  rays  comprising  the  Lg.  Pn  is 
thus  more  vulnerable  to  the  effects  of  non-isotropic  R(0i).  (3)  The  geometrical  spreading 
of  Pn,  even  in  simplified  ID  media,  is  less  stable  than  that  of  the  Lg  (Section  3  of  this 
report;  also  c.f  Sereno  &  Given,  1990,  Zhu  et  al.,  1991). 

§  2.2  Modifications  of  the  algorithm 
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To  reduce  these  bias  for  simultaneous  determination  of  source  Mq,  f^,  and  path  Qq, 
7],  we  have  recently  made  two  modifications  to  the  non-linear  inverse  algorithm  of  Xie 
(1993).  To  present  the  modification,  we  recall  that  in  the  method  of  Xie  (1993),  a  range  of 
combinations  of  Mq  and  values  as  the  possible  solution  are  exhaustively  examined.  For 
each  such  combination,  we  calculate  the  synthetic  source  spectra  5(/)  using  Eq  (2),  and 
thereby  calculating  a  quantity  A''y,  defined  as 


A".  - 
A  ij  - 


Ajifj) 


where  fj  is  the  jth  discrete  frequency.  The  mathematical  expectation  of  ln(A',y)  is 


=  (l-;7,)ln/y-ln(0oi)  , 


(6) 


therefore  we  may  estimate  the  Qi  and  77,  values,  corresponding  to  the  given  Mq  and 
values,  by  a  linear  regression  fitting  using  Eq  (6).  This  enable  us  to  find  one  set  of  possi¬ 
ble  model  vector, 

^  =  (Mq,  fc,Qoi,  TJi,  Qq2,  712,- •  •  Qon,  VnV  »  (7) 

where  N  is  the  total  number  of  stations  recording  the  event.  For  each  possible  set  of  solu¬ 
tions,  we  may  estimate  a  residual,  which  is  defined  as 

^  i=N  j=J(i)  ^ 

Res^  =  E  E  4  ,  (8) 

/=!  j=l 

where  J{i)  is  the  total  number  of  spectral  estimates  available  for  the  ith  station,  and  eij  is 
defined  as 


Eij  =  In  G-HA,)AK/y)  -  In  5(/y)  exp  - 


The  residual,  Res^  is  a  function  of  m.  By  comparing  Res^  calculated  for  all  of  the  possi¬ 
ble  m,  we  find  the  optimal  m  which  gives  the  minimum  of  all  residuals. 

At  the  Ith  station,  the  combined  effect  of  the  station  variable  R{6i),  Xi(f),  Aq,  m  and 
the  instrument  malfiinction/miscalibration  may  raise  or  deduce  A,(/)  by  an  amount  that 
is  most  likely  within  a  factor  of  2  (Xie  &  Mitchell,  1990;  Cong  et  al,  1996).  In  addition, 
the  randomness  in  A(/)  causes  it  to  fluctuate  up  and  down  at  localized  frequencies.  At 
close-in  stations  (say  A;  <  VgQi(f)l(;c /)),  where  the  effect  of  amplitude  reduction  due  to 
finite  Q  is  not  yet  very  pronounced.  A'/,  may  well  be  less  than  zero  at  localized  frequen¬ 
cies  (this  is  the  well-known  "negative  Q"  problem  in  Q  studies).  In  the  method  of  Xie 
(1993),  even  if  negative  A'y  occurs  at  very  localized  frequencies,  mathematical  singulari¬ 
ties  will  occur  when  the  right-hand  side  of  Eq  (5)  is  estimated,  and  the  corresponding  m 
will  then  be  rejected  as  a  reasonable  solution.  This  could  occur  even  when  the  overall 
trend  of  A"(/)  can  still  be  fit  by  a  reasonable  combination  of  Qoi  and  t/,-,  via 
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(10) 


fi-n. 


This  motivated  me  to  make  the  first  modification  of  the  method  of  Xie  (1993).  In  the 
modified  algorithm,  for  a  possible  combination  of  Mq  and  values,  A''y  are  calculated  as 
before.  When  at  any  /th  station  one  or  more  negative  A''y  values  occur,  the  (2oi  and  7, 
values  are  fitted  directly  using  A''y  (Eq  (10)),  rather  than  using  ln(A''y)  (Eq  (6)),  using  a 
iterative,  non-linear  least  square  algorithm  (Eq  (12.159)  of  Aki  &  Richards,  1980).  If  the 
resulting  Qqi  and  7,  values  are  reasonable,  the  corresponding  m  is  saved  as  a  possible 
solution.  In  using  both  synthetic  and  real  Lg  or  Pn  data,  I  have  found  that  when  there  are 
stations  which  are  close  to  the  source  and  are  associated  with  higher  Q  values,  such  that 
Aj  <  VgQi{f)/(7tf),  the  modified  algorithm  is  able  to  find  the  optimal  m  models  that  were 
rejected  by  the  unmodified  method  of  Xie  (1993)  due  to  localized  fluctuations  in  A'y. 

The  second  modification  to  the  algorithm  of  Xie  (1993)  is  that  when  additional  infor¬ 
mation  on  source  spectra  obtained  by  using  the  empirical  Green’s  function  (EGF) 
method,  or  on  two-station  Q  by  measuring  two-station  spectral  ratios  are  available,  the 
spectral  inversion  for  optimal  m  incorporates  these  information.  For  example,  if  an  seis¬ 
mic  event  is  recorded  by  stations  at  varying  distances  such  that  two-station  pairs  may  be 
formed  to  measure  inter-station  Q  values  using  the  standard  spectral  ratio  technique,  then 
the  average  two-station  Q  can  be  obtained  separately,  prior  to  the  simultaneous  non-linear 
spectral  inversion.  In  the  latter  inversion,  a  model  vector  m  that  minimizes  the  residual, 
Res^,  but  fails  to  match  either  the  /<.  values  from  the  EGF  analysis,  or  the  averaged  two 
station  Qq  and  7  values  obtained  earlier,  is  not  considered  to  be  the  optimal  solution.  This 
is  because  that  the  smaller  residual  may  have  been  an  artifact  due  to  non-isotropic  source 
radiation  patterns,  or  due  to  deviations  of  source  spectra  from  the  idealized  stochastic 
models  in  Eq  (2). 

In  an  application  of  the  modified  algorithm  to  the  Lg  excitation  and  propagation 
from  the  1995  western  Texas,  My,  =  5.7  event,  Xie  (1996)  finds  that  these  modifications 
are  very  efficient  in  seeking  for  a  optimal  model  parameter,  m,  that  better  describes  the 
spectral  characteristics  of  the  Lg  excitation  and  propagation  to  over  twenty  broad-band 
stations  in  the  continental  U.S.  I  expect  that  the  modified  algorithm  will  perform  even 
better  for  the  study  of  the  Pn  phase.  This  is  because  (1)  Pn  is  more  affected  by  the  source 
radiation  pattern  due  to  its  simpler  ray  composition.  (2)  Pn  is  more  subjected  to  the  effect 
of  spectral  fluctuations  which  cause  negative  A"y  since  for  Pn,  A,-  is  typically  small  (less 
than  about  12  °)  and  Vg  typically  large  (about  twice  as  that  of  Lg).  These  make  the  con¬ 
dition,  A,  <  V gQi{f)l{7t f),  occur  more  frequently  for  Pn.  (3)  The  correction  of  the  geo¬ 
metrical  spreading  of  Pn  is  much  less  precise  than  that  of  the  Lg.  When  this  correction  is 
overestimated,,  A"^-  will  be  underestimated  and  will  tend  to  be  negative. 
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3.  SYNTHETIC  STUDY  OF  Pn  GEOMETRICAL  SPREADING  IN  CENTRAL 
ASIA 


One  of  the  most  uncertain  factors  in  determining  the  Q  and  source  spectral  parame¬ 
ters  using  Pn  is  its  geometrical  spreading  term  (GST)  (e.g.,  Sereno  et  al,  1988;  Sereno  & 
Given  1990;  Zhu  et  al.,  1991).  Sereno  et  al.  and  Zhu  et  al.  suggested  that  the  best  GST 
for  Scandinavia  and  Canada  are  approximately  described  by  Eq  (3),  with  Aq  and  m  being 
approximately  1  km  and  1.3,  respectively.  We  have  estimated  the  GST  with  synthetic  Pn 
seismograms  for  central  Asia,  using  a  frequency-wave  number  integration  algorithm 
(Saikia,  1994)  and  three  recently  published  velocity  models  for  the  area.  These  models 
are  plotted  in  Figure  1,  they  include:  (1)  model  Ml  of  Roecker  et  al.  (1993)  for  northern 
Tien  Shan,  (2)  the  model  of  Gao  &  Richards  (1994)  for  the  Lop  Nor  test  site,  and  (3)  the 
model  for  Kazakhstan  by  Quin  &  Thurber  (1992).  In  the  calculations,  we  applied  the 
earth  flattening  correction  (Sereno  &  Given,  1988).  Figure  2  shows  the  decay  of  syn¬ 
thetic  Pn  amplitude  spectra  between  0.2  and  2.0  hz,  calculated  using  an  explosive  source 
embedded  at  a  depth  of  1  km,  and  three  velocity  models,  including  (1)  model  Ml  of 
Roecker  et  al.  (1993)  for  northern  Tien  Shan,  (2)  Model  of  Gao  &  Richards  (1994)  for 
the  Lop  Nor  test  site,  and  (3)  Model  for  Kazakhstan  by  Quin  &  Thurber  (1992).  These 
models  are  plotted  in  Figure  1.  In  the  calculations,  we  applied  the  earth  flattening  correc¬ 
tion  as  suggested  by  Sereno  &  Given  (1990).  The  linear  regression  fitting  of  the  loga¬ 
rithm  of  the  synthetic  amplitude  decay  results  in  estimates  of  Aq  and  m  that  are  summa¬ 
rized  in  Table  1.  For  all  of  the  models  used,  the  estimated  Rq  is  subjected  to  large  uncer¬ 
tainties  of  factors  of  25  to  277.  The  uncertainty  in  m  estimates  are  moderate,  ranging 
between  about  0.2  and  0.3.  Taking  these  uncertainties  into  account,  the  Rq  and  m  values 
in  Table  1  agree  with  the  values  by  Sereno  et  al.  (1988)  and  Zhu  et  al.  (1991)  for  Scande- 
navia  and  Canada.  In  Figure  2,  we  also  indicate  the  local  deviations  of  the  synthetic 
amplitude  from  the  regressional  fit.  For  all  of  the  three  model  used,  one  should  expect  a 
local  variation  of  about  a  factor  of  2.  This  is  in  contrast  to  what  we  found  for  synthetic  Lg 
amplitude  decay  (not  shown  in  this  report),  in  which  the  local  deviation  is  only  in  the 
order  of  10%.  Therefore,  the  large  uncertainties  in  Rq,  m  estimates  and  the  large  fluctua¬ 
tion  in  Pn  amplitude  decay  from  the  regressional  trend  clearly  indicate  that  an  attempt  to 
precisely  model  Pn  GST,  even  in  the  simplified  ID  cases,  is  unsuccessful.  In  practice, 
this  attempt  is  even  more  unsuccessful  if  we  take  into  account  that  the  real  structure 
varies  three  dimensionally.  As  of  this  writing.  We  plan  to  used  an  m  of  1.3  and  Rq  of  1 
km  (Sereno  et  al.,  1988)  throughout  this  research. 
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Table  1.  Values  of  Aq  and  m  from  Pn  Synthetics 


Velocity  Model  Used 

Ao 

m 

Ml  (Roeckereta/.) 

23.42  (x/29)km 

1.43  ±0.24 

Gao  &  Richards  (1994)* 

24.1  (x/25)km 

1.37  ±0. 17 

Quin  &  Thurber 

14.1  (X/  277)  km 

1.36  ±0.29 

*  Used  model  Ml  for  mantle  below  75  km. 
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Vp  models 
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et  al.  (1993).  The  dashed  line  is  the  model  of  Gao  &  Richards  (1994)  for  the  Lop  Nor  Test  Site,  combined  with  model  Ml  at  depths 
>  75  km.  The  dotted  line  is  the  model  by  Quin  &  Thurber  (1992)  for  Kazakhstan. 


4.  Pn  AMPLITUDE  VARIATION  L\  CENTRAL  ASIA 

We  have  collected  more  than  one  hundred  Pn  spectra  recorded  at  16  broad-band 
IRIS,  Kyrghistan  network  (KNET),  Kazakhstan  Network  (KZN)  and  CDSN  stations, 
from  three  underground  nuclear  explosions  and  14  earthquakes  in  central  Asia.  A  wealth 
of  information  have  been  gained  on  the  characteristics  of  excitation  and  propagation  of 
Pn  in  central  Asia.  Quantitative  inversion  for  Mq,  fc  and  path  Q  have  not  been  conducted 
to  most  of  the  spectra  obtained,  since  we  spent  most  of  the  last  year  for  algorithm  and 
software  development,  as  well  as  for  understanding  the  Pn  geometrical  spreading.  We 
have,  however,  extensively  examined  feamres  of  time  domain  and  frequency  domain  Pn 
amplitude  from  many  earthquakes  and  several  explosions.  One  important  finding  if  that 
the  Pn  amplitude  from  Lop  Nor  explosions  vary  drastically  across  the  KNET  stations.  As 
an  example.  Figure  3  shows  the  record  section  from  the  Oct.  7,  1994  Lop  Nor  explosion, 
which  is  located  to  the  east  of  the  KNET  stations  (Figure  5).  The  Pn  amplitude  varies  by 
a  factor  of  about  19  across  the  KNET,  despite  the  relative  small  span  (~  2°)  of  the  net¬ 
work.  Figure  4  shows  the  corresponding  amplitude  spectra  of  Pn,  obtained  using  a  win¬ 
dow  with  an  effective  length  of  4.5  s.  The  amplitude  variation  primarily  occur  at  low  (- 
Ihz)  frequencies.  For  instance,  the  1  hz  amplitude  is  a  factor  of  27  higher  at  station  KBK 
than  at  AML. 

Various  attempts  have  been  made  to  explore  the  origin  of  this  amplitude  variation. 
We  studied  frequency-dependent  polarization  of  Pn  across  KNET  from  Lop  Nor  explo¬ 
sions,  but  found  that  the  polarization  is  strongly  affected  by  the  topography  near  individ¬ 
ual  stations,  thus  is  not  very  useful  in  inferring  the  origin  of  the  amplitude  variation.  We 
studied  Pn  amplitude  across  KNET  from  earthquakes  in  southern  Tien  Shan,  and  found 
that  the  pattern  of  Pn  amplitude  variation  differs  from  that  of  Lop  Nor  explosions.  Figure 
6  shows  a  record  section  from  a  southern  Hen  Shan  earthquake.  In  Figure  6,  Pn  ampli¬ 
tudes  at  stations  KBK  and  AML  are  low  and  high,  respectively,  opposite  to  the  trend 
observed  in  Figure  3.  Therefore,  a  systematic  site  response  is  ruled  out  as  the  cause  of  the 
amplitude  variation.  We  also  analyzed  Pn  from  earthquakes  in  northern  Xinjiang,  west  of 
Lake  Bosten  Hu  (Figure  5).  Figure  7  shows  the  record  section  for  one  of  such  earth¬ 
quakes.  Generally,  Pn  is  very  weak  across  the  KNET,  but  the  pattern  of  amplitude  varia¬ 
tion  is  similar  to  that  found  from  Lop  Nor  explosions:  Pn  amplitude  reaches  its  maximum 
at  stations  KBK  and  KZA  and  drops  to  a  minimum  at  more  distant  stations.  This  suggests 
that  some  3D  velocity  anomaly  along  the  east-west  oriented  paths  is  responsible  for  the 
amplitude  variation.  Figure  8  shows  a  synthetic  Pn  ray  pattern  through  an  ID  model, 
which  is  adapted  from  model  Ml  (last  section)  by  adding  a  velocity  gradient  of  2  x  10~^ 
s/km  (before  the  earth  flattening  correction)  into  the  uppermost  mantle  between  depths  of 
50  to  150  km.  This  ray  pattern  is  stable  and  unrealistic,  since  the  observed  large  ampli¬ 
tude  variation  must  be  associated  with  some  kind  of  unstable  ray  pattern.  The  easiest  way 
to  introduce  an  unstable  pattern  is  to  introduce  a  3D  topography  to  the  Moho,  across 
which  the  velocity  contract  is  the  most  drastic  for  all  boundaries  at  depth.  Figure  9  shows 
ray  pattern  in  a  2D  strucmre,  which  is  adapted  from  the  ID  structure  in  Figure  8  by 
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the  Pn  phase.  The  time  corresponding  to  a  group  velocity  of  3.5  km/s  (typical  of  Lg)  is  marked  by  a  straight  line.  The  traces  are 
normalized  by  a  common  factor.  Station  name  and  peak  Pn  and  Lg  amplitudes  are  indicated  near  each  trace.  Note  that  at  KBK,  the 
amplitude  is  12  times  of  that  at  a  nearby  station  (CHM),  and  19  times  of  that  at  station  AML. 


MAK,  840,0  km  ULHL,  1060.6  km  KZA,  1149.1  km  KBK,  1173.7  km 


apninduiy  apn^iiduiy 
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Figure  4.  Pn  and  noise  spectra  from  the  Oct.  7,  1994  Lop  Nor  explosion,  calculated  for  ten  KNET  and  Kazakhstan  stations  (all  but 
the  lower  right  panels),  with  an  effective  window  length  of  4.5  s.  The  lower  right  panel  is  the  spectral  ratio  of  KBK/AML.  Note  the 
wide  range  of  the  amplitudes  at  1  Hz  for  various  stations,  and  the  factor  of  27  at  1  Hz  in  the  KBK/AML  spectral  ratio. 
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representative  KNET  station  (station  KBK)  are  plotted. 


KNET  Record  Sectio 
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phase.  Plotted  in  the  same  manner  as  in  Figure  2.  Note  that  the  Pn  amplitude  at  KBK  is  lower  than  at  AML. 
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phase.  Plotted  in  the  same  manner  as  in  Figures  2  and  5. 


purpose  since  the  velocity  gradient  used  is  likely  larger  than  the  average  value  in  central  Asia. 
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tance,  the  ray  density  shows  a  pattern  of  normal-high-low;  this  corresponds  to  the  observed  amplitude  variation  in  Figures.  3,  4  and 


introducing  a  2D  topography  to  Moho.  The  depth  of  Moho  in  Figure  9  varies  from  55  km 
to  45  km  between  distances  of  about  900  and  1000  km.  This  simple  Moho  topography 
results  in  an  amplitude  pattern  that  is  qualitatively  in  agreement  with  that  observed  for  the 
explosions  and  earthquakes  from  the  east  of  the  KNET  (Figures  3  and  7):  the  Pn  ray  den¬ 
sity  (which  is  proportional  to  the  amplitude)  remains  normal  until  the  distance  reaches 
about  960  km,  or  some  60  km  away  from  the  point  where  the  Moho  thinning  begins.  Pn 
ray  density  at  around  the  distance  of  960  km  increases  drastically.  As  distances  further 
increases,  the  Pn  ray  density  drops  to  abnormally  low.  This  example  demonstrates  how  a 
2D  or  3D  Moho  topography  can  drastically  change  the  Pn  amplitude.  There  may  be  other 
3D  velocity  structures  that  could  produce  similar  Pn  amplitude  variations.  Any  3D  Moho 
topography  responsible  for  the  drastic  amplitude  variation  may  well  be  associated  with 
the  drastic  variations  in  surface  topography  and  geology  from  Lake  Issyk-kul  to  the 
mountains  in  the  west  (Figure  5).  Unfortunately,  the  Pn  data  that  we  retrieved  so  far  is  not 
sufficient  to  uniquely  constrain  the  real  3D  structure. 


5.  PRELIMmARY  EVALUATION  OF  Pn/Lg  RATIO  AS  A  DISCRIMINANT 

Comparing  Figures  3  and  7,  one  would  be  tempted  to  draw  an  easy  conclusion  that 
the  Pn/Lg  amplitude  ratio  works  well  as  a  discriminant  between  explosions  and  earth¬ 
quakes.  However,  this  conclusion  is  premature.  The  reason  for  the  low  Pn  amplitudes  in 
Figure  7  is  likely  caused  by  a  null  in  the  source-to-Pn  radiation  pattern  since  the  earth¬ 
quakes  occurred  east  of  Lake  Bosten  Hu,  where  earthquakes  sources  are  typically  charac¬ 
terized  by  an  east- west  strike-slip  motion  (Center  for  Analysis  and  Prediction,  1991).  The 
earthquakes  from  southern  Tien  Shan  generate  stronger  Pn  signals  across  KNET  (Figure 
6),  suggesting  that  the  Pn/Lg  ratio  from  earthquakes  is  significantly  affected  by  the  focal 
mechanism.  Therefore,  any  attempt  to  use  Pn/Lg  amplitude  ratio  as  a  discriminant  must 
account  for  the  instabilities  of  the  Pn  amplitude  (and  to  a  less  extent,  Lg  amplitude) 
caused  by  3D  velocity  complications,  as  well  as  the  source-to-Pn  radiation  pattern.  We 
are  currently  analyzing  Pn  from  more  earthquakes  to  quantify  these  instabilities. 
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6.  MEASUREMENT  OF  Mq,  A  AND  PATH  Q  USING  Pn  FROM  THE  OCTOBER 
7, 1994  LOP  NOR  EXPLOSION 

In  a  preliminary  analysis,  we  measured  the  source  and  path  spectral  parameters  using 
Pn  from  the  Oct.  7,  1994  Lop  Nor  explosion,  with  the  modified  algorithm  mentioned  in  a 
previous  section,  and  the  modified  Mueller-Murphy  source  model  (Sereno  et  al.,  1988), 
with  an  overshoot  parameter  {P  in  Eq  (2))  of  1.0.  Station  KBK  is  not  used  due  to  the 
abnormal  large  amplitude  (Figures  3,4  and  7).  The  resulting  Mq  and  values  are  3.5 
X  10^^  dyne-cm  and  1.6  Hz,  respectively.  For  comparison,  a  previous  inversion  using  Lg 
(Xie  et  al,  1996)  has  resulted  in  an  Mq  of  5.5  x  10^^  dyne-cm,  and  an  of  0.8  Hz.  The 
averaged  Pn  Qq  and  77  values  are  379  and  0.7,  respectively,  for  the  nine  paths  connecting 
the  Lop  Nor  and  stations  MAK,  ULHL,  KZA,  CHM,  UCH,  USP,  EKS,  AML,  KUR.  All 
but  the  first  of  these  stations  belong  to  the  KNET.  The  lowest  and  highest  Qq  values  are 
204  (to  station  AML)  and  573  (to  station  KZA).  Figure  10  shows  the  synthetic  source 
spectra  for  that  event,  versus  the  observed  spectra  after  path  Q  correction.  The  overall  fit 
between  the  synthetic  source  spectrum  and  the  station  average  of  the  observed  (lower 
right  panel)  is  rather  good. 


7.  SOFTWARE  DEVELOPMENT 

The  computer  softwares  developed  during  the  last  year  include: 

1.  Program  "Pn96"  which  takes  the  time  series  containing  Pn  or  any  seismic  phases, 
windows  the  time  series  between  any  two  specified  lapse  times  or  group  velocities 
(adjustable),  and  calculate  Fourier  spectra  of  the  windowed  phase  and  noise.  It  then  cor¬ 
rect  the  effect  of  an  specified  geometric  spreading  and  reduce  the  spectra  to  source.  This 
program  is  written  for  studying  Pn,  Pn  coda  and  Pg,  but  can  be  used  for  studying  any 
regional  phases.  Figure  1 1  shows  and  example  of  the  application  of  this  program. 

2.  Program  "PnQ96"  which  takes  the  spectral  output  of  Pn96.f  from  multiple  sta¬ 
tions,  and  searches  though  all  of  the  possible  combinations  of  Mq,  fc  (and  thereby  m  in 
Eq  (7))  to  find  the  optimal  m.  The  Qq  and  77  are  allowed  to  be  station-variable.  Alterna¬ 
tively,  Qo  7]  can  be  assumed  to  be  station  invariant  in  which  case  site  response,  X(f), 
is  assumed  to  be  station- variable.  If  there  is  any  information  on  the  source  spectrum  from 
the  EGF  method,  or  on  path  Q  from  two-station  method,  this  program  incorporate  these 
information  during  the  search.  During  the  search,  if  negative  A'-  in  Eq  (5)  occurs,  this 
program  finds  Qqi  and  77,  by  an  interactive,  non-linear  algorithm. 

3.  A  set  of  programs  that  finds  the  average  two-station  Q  using  spectral  output  of 
"Pn96"  by  stacking  spectral  ratios  for  all  pairs  of  two  stations  available. 
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Figure  10.  Synthetic  Pn  source  spectra  for  KNET  and  IRIS  stations  recording  the  Oct.  7,  1994,  Lop  Nor  explosion,  versus  the 
observed.  The  lower  right  panel  is  the  station  average.  The  synthetic  spectra  are  calculated  using  a  yd  of  1.0  and  the  optimal  source 
spectral  parameters  (Mq  =  3.5  x  10^^  dyne-cm,  4=1.6  Hz)  obtained  in  the  inversion.  The  distance,  Qq,  t?  and  values  are  written  on 
the  top  of  the  panels.  The  unit  of  the  amplitude  is  the  same  as  that  of  Figure  3  of  Xie  et  al  (1996). 
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Figure  11.  An  example  of  the  application  of  program  Pn96.  The  top  panel  is  the  time  series  con- 
taining  Pn  from  a  Lop  Nor  explosion  at  a  KNET  station.  Group  velocities  and  time  after  the 
beginning  of  the  window  are  indicated  on  the  top  and  bottom  of  the  panel.  The  smooth  curve 
represents  a  20  percent  cosine  window  used  in  the  analysis.  The  bottom  panel  show  the  Fourier 
amplitude  spectra  of  Pn  and  noise. 
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8.  CONCLUSIONS  AND  FUTURE  WORKS 


Much  of  the  first  year  of  this  research  has  been  devoted  to  software  development, 
data  collection  and  processing,  and  synthetic  studies  of  the  Pn  propagation.  Several 
important  findings  have  been  made  on  the  characteristics  of  Pn  excitation  and  propagation 
in  central  Asia.  These  include: 

(1)  Two  modifications  to  the  spectral  inverse  algorithm  of  Xie  (1993)  serves  well  to 
improve  our  ability  to  seek  for  more  reliable  and  robust  estimates  of  source  Mq, 
and  path  Qq,  tj  values  using  regional  wave  spectra. 

(2)  Pn  geometrical  spreading  estimated  using  synthetics  in  ID  structures  for  central  Asia 
suggests  that  it  agrees  roughly  with  a  trend  of  although  the  localized  fluctua¬ 
tion  is  quite  significant  (by  a  factor  of  2  or  more)  from  that  trend. 

(3)  The  observed  Pn  amplitude  varies  drastically  (by  a  factor  of  up  to  20  to  30)  across 
KNET  from  Lop  Nor  explosions  and  Northern  Xinjiang  earthquakes,  which  are 
located  to  the  east  of  the  stations.  This  amplitude  variation  primarily  occurs  at  lower 
(~  1  hz)  frequencies,  and  does  not  correlate  with  the  amplitude  variation  across 
KNET  from  events  in  the  south.  The  most  likely  cause  of  the  amplitude  variation  is  a 
short  wavelength  variation  in  the  Moho  topography. 

(4)  Pn/Lg  amplitude  ratio  also  appears  to  be  highly  dependent  on  the  source-to-Pn  radia¬ 
tion  pattern.  This  ratio  can  be  used  efficiently  to  identify  earthquakes  when  the  sta¬ 
tion  is  located  at  the  Pn  radiation  null.  For  other  source-station  configurations,  how¬ 
ever,  the  efficiency  of  the  ratio  as  a  discriminant  between  explosions  and  earthquakes 
needs  further  investigation. 

(4)  Preliminary  inversion  of  a  well-recorded  Lop  Nor  explosion  suggests  that  the  Mq 
and  fc  values  derived  using  Pn  are  both  higher  than  those  derived  using  Lg.  Inver¬ 
sion  of  Pn  spectra  from  many  more  events,  with  an  improved  method,  is  underway. 

Research  conducted  so  far  suggests  that  any  attempts  to  use  Pn/Lg  amplitude  ratio  as 
a  discriminant,  or  to  invert  for  source/path  spectral  parameters  must  take  into  account  the 
Pn  amplitude  variations  caused  by  effects  of  3D  velocity  complications  and  the  source 
radiation  pattern.  Since  both  effects  are  supposed  to  be  reduced  when  Pn  coda  is  used,  we 
have  already  implemented  an  algorithm  to  invert  for  Mq,  path  Qq  and  77  using  the 
entire  wavetrain  that  includes  Pn,  Pn  coda  and/or  to  Pg.  We  will  compare  the  stability  of 
the  results  using  Pn  alone  with  those  using  the  averaged  spectra  of  Pn  and  Pn  coda.  Since 
almost  all  necessary  software  development  has  been  achieved,  the  next  phase  of  the 
research  will  be  mainly  devoted  to  data  processing,  inversion  and  interpretation. 

By  collecting  Pn  records  at  KNET  stations  from  more  Northern  Xinjiang  earth¬ 
quakes,  particularly  those  earthquakes  whose  focal  mechanisms  differ  from  the  typical 
east-west  strike  slip,  we  will  continue  to  explore  the  cause  for  the  large  Pn  amplitude 
variation  across  the  KNET  for  the  east-west  oriented  paths.  In  particular,  the  existence  of 
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any  3D  Moho  topography  and  its  link  to  the  topographic/geological  variation  from  Lake 
Issyk-kul  to  Northern  Tien  Shan  will  be  investigated. 
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PACIFIC  NORTHWEST  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  999,  MS  K6-40 
RICHLAND,  WA  99352 


PACIFIC  NORTHWEST  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  999,  MS  K7-22 
RICHLAND,  WA  99352 


PACIFIC  NORTHWEST  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  999,  MS  K5-72 
RICHLAND,  WA  99352 


PACIFIC  NORTHWEST  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  999,  MS  K6-84 
RICHLAND,  WA  99352 


PACIFIC  NORTHWEST  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

PO  BOX  999,  MS  K5-12 
RICHLAND,  WA  99352 


FRANK  PILOTTE 
HQ/AFTAC/TT 
1030  S.  HIGHWAY  AI A 
PATRICK  AFB,  FL  32925-3002 


KEITH  PRIESTLEY 
DEPARTMENT  OF  EARTH  SCIENCES 
UNIVERSITY  OF  CAMBRIDGE 
MADINGLEY  RISE,  MADINGLEY  ROAD 
CAMBRIDGE,  CB3  OEZ  UK 


JAY  PULLI 

RADIX  SYSTEMS,  INC. 
6  TAFT  COURT 
ROCKVILLE,  MD  20850 


PAUL  RICHARDS 
COLUMBIA  UNIVERSITY 
LAMONT-DOHERTY  EARTH  OBSERVATORY 
PALISADES,  NY  10964 


DAVID  RUSSELL 
HQ  AFTAC/TTR 
1030  SOUTH  HIGHWAY  AI  A 
PATRICK  AFB,  FL  32925-3002 


CHANDAN  SAIKIA 

WOOODWARD-CLYDE  FEDERAL  SERVICES 
566  EL  DORADO  ST.,  SUITE  100 
PASADENA,  CA  91 101-2560 


SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  61 16 

MS  0750,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0750 

SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  9311 

MS  1 159,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-1 159 

SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 
DEPT.  5736 

MS  0655,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0655 

THOMAS  SERENO  JR. 

SCIENCE  APPLICATIONS  INTERNATIONAL 

CORPORATION 

10260  CAMPUS  POINT  DRIVE 

SAN  DIEGO,  C A  92121 


SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

DEPT.  5704 

MS  0979,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0979 

SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

DEPT.  5791 

MS  0567,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0567 

SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

DEPT.  5704 

MS  0655,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0655 

SANDIA  NATIONAL  LABORATORY 
ATTN:  TECHNICAL  STAFF  (PLS  ROUTE) 

DEPT.  6116 

MS  0750,  PO  BOX  5800 
ALBUQUERQUE,  NM  87185-0750 

AVI  SHAPIRA 
SEISMOLOGY  DIVISION 

THE  INSTITUTE  FOR  PETROLEUM  RESEARCH  AND 
GEOPHYSICS 

P.O.B.  2286,  NOLON  58122  ISRAEL 
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ROBERT  SHUMWAY 
410MRAK  HALL 
DIVISION  OF  STATISTICS 
UNIVERSITY  OF  CALIFORNIA 
DAVIS,  CA  95616-8671 

DAVID  SIMPSON 
IRIS 

1616  N.  FORT  MEYER  DRIVE 
SUITE  1050 

ARLINGTON,  VA  22209 

BRIAN  SULLIVAN 
BOSTON  COLLEGE 
INSITUTE  FOR  SPACE  RESEARCH 
140  COMMONWEALTH  AVENUE 
CHESTNUT  HILL,  MA  02167 

NAFI TOKSOZ 

EARTH  RESOURCES  LABORATORY,  M.I.T. 
42  CARLTON  STREET,  E34-440 
CAMBRIDGE,  MA  02142 


GREG  VAN  DER  VINK 
IRIS 

1616  N.  FORT  MEYER  DRIVE 
SUITE  1050 

ARLINGTON,  VA  22209 

TERRY  WALLACE 
UNIVERSITY  OF  ARIZONA 
DEPARTMENT  OF  GEOSCIENCES 
BUILDING  #77 
TUCSON,  AZ  85721 

JAMES  WHITCOMB 
NSF 

NSF/ISC  OPERATIONS/EAR-785 
4201  WILSON  BLVD.,  ROOM785 
ARLINGTON,  VA  22230 

JIAKANG  XIE 
COLUMBIA  UNIVERSITY 
LAMONT  DOHERTY  EARTH  OBSERVATORY 
ROUTE  9W 

PALISADES,  NY  10964 

OFFICE  OF  THE  SECRETARY  OF  DEFENSE 
DDR&E 

WASHINGTON,  DC  20330 


TACTEC 

BATTELLE  MEMORIAL  INSTITUTE 
505  KING  AVENUE 

COLUMBUS,  OH  43201  (FINAL  REPORT) 


MATTHEW  SIBOL 
ENSCO,  INC. 

445  PINEDA  COURT 
MELBOURNE,  FL  32940 


JEFFRY  STEVENS 
MAXWELL  TECHNOLOGIES 
P.O.  BOX  23558 
SAN  DIEGO,  CA  92123 


DAVID  THOMAS 
ISEE 

29100  AURORA  ROAD 
CLEVELAND,  OH  44139 


LAWRENCE  TURNBULL 

ACIS 

DCI/ACIS 

WASHINGTON,  DC  20505 


FRANK  VERNON 

UNIVERSITY  OF  CALIFORNIA,  SAN  DIEGO 
SCRIPTS  INSTITUTION  OF  OCEANOGRAPHY  IGPP,  0225 
9500  GILMAN  DRIVE 
LA  JOLLA,  CA  92093-0225 

DANIEL  WEILL 
NSF 

EAR-785 

4201  WILSON  BLVD.,  ROOM  785 
ARLINGTON,  VA  22230 

RU  SHAN  WU 

UNIVERSITY  OF  CALIFORNIA  SANTA  CRUZ 
EARTH  SCIENCES  DEPT. 

1 156  HIGH  STREET 
SANTA  CRUZ,  CA  95064 

JAMES  E.  ZOLLWEG 
BOISE  STATE  UNIVERSITY 
GEOSCIENCES  DEPT. 

1910  UNIVERSITY  DRIVE 
BOISE,  ID  83725 

DEFENSE  TECHNICAL  INFORMATION  CENTER 
8725  JOHN  J.  KINGMAN  ROAD 
FT  BELVOIR,  VA  22060-6218  (2  COPIES) 


PHILLIPS  LABORATORY 
ATTN:  XPG 
29  RANDOLPH  ROAD 
HANSCOM  AFB,  MA  01731-3010 
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PHILLIPS  LABORATORY 
ATTN:  GPE 
29  RANDOLPH  ROAD 
HANSCOM  AFB,  MA  0173 1-3010 


PHILLIPS  LABORATORY 
ATTN:  TSML 
5  WRIGHT  STREET 
HANSCOM  AFB,  MA  01731-3004 


PHILLIPS  LABORATORY 
ATTN:  PL/SUL 
3550  ABERDEEN  AVE  SE 
KIRTLAND,  NM  871 17-5776  (2  COPIES) 
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